Approximating invariant densities of metastable 

systems 

Cecilia Gonzalez Tokman* Brian R. Hunt* Paul Wright* 

May 2, 2009 

Abstract 

We consider a piecewise smooth expanding map of the interval pos- 
sessing two invariant subsets of positive Lebesgue measure and ex- 
actly two ergodic absolutely continuous invariant probability measures 
(ACIMs). When this system is perturbed slightly to make the invari- 
ant sets merge, we describe how the unique ACIM of the perturbed 
map can be approximated by a convex combination of the two initial 
ergodic ACIMs. 

1 Introduction 

Metastable systems are studied in relation with phenomena ranging from 
molecular [MDHS06] to oceanic [FPET07] dynamics. Typical trajectories 
of these systems remain in one of its almost invariant (metastable or quasi- 
stationary) components for a relatively long period of time, but eventually 
switch to a different component and repeat this behavior. Quantitative as- 
pects of these phenomena have been studied through eigenvalue and eigen- 
vector approximation techniques for Markov models [MSF05, FP08]. Here, 
we are concerned with rigorous approximation results for eigenvectors-in 
particular those that correspond to stationary measures of the dynamics-in 
a more general (non-Markov) setting. 

Broadly, our setting concerns the approximation of absolutely continuous 
invariant probability measures (ACIMs) for certain hyperbolic maps with 
metastable states. These systems arise from perturbing an initial system 
To with two disjoint invariant sets //, I r of positive Lebesgue measure. The 
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initial map has two mutually singular ergodic ACIMs, m and [i r . When 
To is perturbed in such a way that i/ and I r lose their invariance and the 
perturbed map T t has only one ACIM fj, e , we are interested in approximating 
fi e using m and \i r . Specifically, the systems we consider are piecewise C 2 
expanding maps of an interval; see Figure 1. 

Our results can be understood in the context of dynamical systems with 
holes as follows. As the invariance of the two initially invariant sets is 
destroyed by the perturbation, we think of the small set of points // n T~ I r 
that switch from I; to I r , and likewise the set I r n T~ I[, as being holes in 
the initially invariant sets. From this point of view we expect to be able to 
approximate fj, e , for small e, by a convex combination ol\i\ + (1 — a)[i r of the 
two initially invariant measures, with the ratio a /{I — a) depending on the 
relative sizes of the holes. 
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Figure 1: Dashed: initial system. Thick: metastable system. 



Before discussing our results, we present two illustrative examples. We 
begin with a simple random system. Consider the family of Markov chains 
in two states I and r, with transition matrices 

where e = (e/^. r , e r ^i). We are interested in the behavior when e ~ 0. When 
e = 0, the two sets I\ = {1} and I r = {r} are invariant, giving rise to the 
two ergodic stationary probability measures /i/ = 5i and /u r = S r . When 
ti-^r > 0, there is a unique stationary probability measure 

He = oteM + (1 - a e W, where — — — = e ^l. 

1 - a e ei-^ r 
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Observe that the ratio of the weights a e /(l — a e ), i.e. ^ e (I;)//x e (I r ), is equal 
to the inverse ratio of the sizes of the holes, e r ^/e^ r . 

Next, we consider two billiard tables T>[,T> r in the plane, as indicated in 
Figure 2. For * £ {/,r}, let T* : I* (3 be the corresponding billiard map, 
i.e. the Poincare map for the first return of the billiard flow to dT>*. We 
use \dT>*\ to denote the perimeter of D*. A general reference for hyperbolic 
billiards is [CM06], where one can find the background for the assertions 
below. We use the usual coordinates (s,ip) on /*, where s is arc length 
on dV*, and cp £ [— tt/2, +vr/2] is the angle between the outgoing velocity 
vector and the inward pointing normal vector to dT>* . Then it is well known 
that T* leaves (normalized) Liouville measure /i* invariant, where /i* has 
the density := d^/dsdip = [2 cos ip. Next, for e > 0, let h e 

be a subsegment of dT>i n dV r of length e, and let V t be the billiard table 
resulting after h t is removed. The corresponding density for the invariant 
Liouville measure of the billiard map is </> e = [2(|<9D;| + \dV r \ — 2e)] _1 cos ip. 
Thus as e — > 0, 

<£ e ->■ «0/ + (1 - a)(/> r , where = Jf?' 

1 — a | aP r | 

provided some care is taken to define all of the density functions involved 
on the same space. Note that if we define the holes H* :t := T~ 1 (h e x 
[— it/2, +7r/2]), then we can rewrite a/(l — a) = fj, r (H rtt )/ /ii(Hi ;e ), so that 
again the ratio of the weights equals the inverse ratio of the sizes of the 
holes. This example is most meaningful when T/, T r , and T e are all ergodic, 
which is the case for the tables in Figure 2; see §8.15 in [CM06]. 




Figure 2: Two ergodic billiard tables connected by a hole. 

In our main result, Theorem 1, we show a corresponding result in the 
deterministic setting of piecewise C 2 expanding maps, under fairly general 
conditions described in §2. We show that as e — > the invariant density 
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(j) t of T e converges in L to a convex combination of the ergodic invariant 
densities of To, with the ratio of the weights given by the limiting inverse 
ratio of the sizes of the holes. We emphasize that our results do not require 
any of the piecewise expanding maps involved to have a Markov partition. 

The density <j) t corresponds to an eigenvector with eigenvalue 1 for the 
Perron-Frobenius operator acting on a suitable space of functions. Our 
assumptions imply that for e > 0, the operator has 1 as a simple eigenvalue, 
and also another real simple eigenvalue slightly less than 1. In Theorem 2, 
we characterize the eigenvectors of this lesser eigenvalue by showing that 
asymptotically they lie on the line spanned by dpi/dx — d\i r jdx. 

Unlike the two examples above, in the setting of piecewise C 2 expanding 
maps we have no explicit formulas for the invariant densities; even their 
existence is nontrivial. Our methods rely on the fact that the densities 
of the ACIMs for T e are of bounded variation [LY73]. Hence, they can 
be decomposed into regular and singular (or saltus) parts, as in [Bal07]. 
The key technical portions of our proofs include estimating and exploiting 
the locations and sizes of the jumps at the discontinuities of the invariant 
densities, which occur on the forward trajectories of the critical points of T e . 

This work is related to other recent work involving metastable systems 
and piecewise expanding maps. Recently, [ ] studied metastable systems 
arising from piecewise smooth uniformly expanding maps with two invari- 
ant intervals. They perturbed such an initial map by a family of Markov 
operators close to the identity to produce a family of metastable systems 
for e > 0. The associated Perron-Frobenius operators acting on a suitable 
space of functions have 1 as a simple eigenvalue and another simple eigen- 
value p e < 1. As e — ► 0, p e — ► 1, and the authors rigorously computed the 
derivative lim e _ +0 +(l — Pe)/ e - This provides information on the stationary 
exchange rate between the metastable states. Their work may be used to 
show a corresponding result in our setting. 

Our work is also related to current and ongoing investigations on linear 
response. These problems have the feature that, as Ruelle [Ruc98] puts it, 
it is possible to formulate conjectures based on intuition or formal calcula- 
tions, but the proofs often involve overcoming intricate technicalities. In our 
setting, we know that fJL e {Ii) — > a as e — ► 0. A pertinent open problem would 
be to try and characterize the higher-order terms 7Z(e) := (J, e (Il) ~ a - We do 
not expect TZ(e) to be differentiable at e = in general. As shown in [BS08], 
linear response fails precisely when the perturbations T e are transverse to 
the topological class of To, at least for certain piecewise expanding unimodal 
maps To that are topologically mixing. Results of [Kel82] show that in that 
setting the unique ACIM 4> e of T e satisfies \4> t — 4>o\ L i = O(eloge), where 
4>o is the unique ACIM of Tq. [Bal07] gives examples where this estimate is 
optimal. 

Another problem for further research is extending our results to higher 
dimensional piecewise hyperbolic maps. While we use techniques specific 
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to one- dimensional maps, we are optimistic that the main elements of our 
proof, found in §3.2, can be generalized. 

2 Statement of results 

In this section, we define a class of dynamical systems with two nearly in- 
variant (metastable) subsets. They are perturbations of a one-dimensional 
piecewise smooth expanding map with exactly two invariant subintervals 1\ 
and I r of positive Lebesgue measure. On each of these intervals, the unper- 
turbed system has a unique ACIM. The perturbations break this invariance 
by introducing what we consider to be holes in the intervals; the hole(s) in Ii 
map to I r and vice versa. Each perturbed system will have only one ACIM, 
and we will determine an asymptotic formula for its density in terms of the 
invariant densities of the unperturbed system. 

Let I = [0, 1]. In this paper, a map T : I O is called a piecewise C 2 
map with C = {0 = Co < c\ < ■ ■ • < = 1} as a critical set if for each 
i, T\r c . jC . ) extends to a C 2 function on a neighborhood of [c^Cj+i]. We 
call T uniformly expanding if its minimum expansion, inf^g/y^ |Tq(:e)|, is 
greater than 1. As is customary for piecewise smooth maps, we consider 
T to be bi- valued at points c, G C where it is discontinuous. In such cases 
we let T(q) be both values obtained as x approaches q from either side, 
and T(cj-t) the corresponding right and left limits. If a, 6 € C, T\t a u will be 
used to specifically denote the restriction of T with T|[ ab j(a) = T(a+) and 
T| M (6)=T(6„). 

We use Leb to denote normalized Lebesgue measure on I and L 1 to 
denote the space of Lebesgue integrable functions on /, with norm \f\ L i = 
Ji \f( x )\ dx. Also, for / : I — > C, we let {fl^ be the supremum of / over / 
and var(/) be the total variation of / over /; that is, 

n 

var(/) = sup{^ \f(xi) - : n > 1, < x < x x < ■ ■ • < x n < 1}. 

i=l 

For clarity of presentation, we do not state our results under the broadest 
possible assumptions. However, see §2.4 for a number of relaxations of the 
hypotheses below. 

2.1 The initial system and its perturbations 

We assume that the unperturbed system is a piecewise C 2 uniformly ex- 
panding map Tq : I O with Co = {0 = co,o < c^o < • • • < c^o = 1} as a 
critical set. There is a boundary point b £ (0, 1) such that // := [0,6] and 
I r := [b, 1] are invariant under To, i.e. for * G {l,r}, 7b|/,(/*) C I*. The ex- 
istence of an ACIM of bounded variation for To | j t is guaranteed by [LY73] . 
We assume in addition: 
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(11) Unique ACIMs on the initially invariant set. 

To|/„ has only one ACIM whose density is denoted by (ft* := d/i*/dx. 

The uniqueness of such an ACIM can be guaranteed by transitivity or by 
additional conditions described in [LY78]. From (II), it follows that all 
ACIMs of To are convex combinations of the ergodic ones, fi[ and ji r . 

We define the points in Hq := T _1 {6} \ {b} to be infinitesimal holes. 
These are all points that map to the boundary point b, except possibly b 
itself. Our reasons for excluding b from the set of infinitesimal holes will 
be explained in §2.4. An immediate consequence of this definition is that 
Hq C Co. 

(12) No return of the critical set to the infinitesimal holes. 
For every k> 0, (T fc C ) n H = 0. 

As we will see in §4.2, this implies that (ft* is continuous at each of the 
infinitesimal holes in I*. 

(13) Positive ACIMs at infinitesimal holes. 

<fti is positive at each of the points in Hq n , and <ft r is positive at each 
of the points in Hq n I r ■ 

For example, this will be the case if To\j t and ?o|j r are weakly covering, 1 
see [Liv95]. 

(14) Restriction on periodic critical points. 
Either 

(I4a) mf xeI \ Co m(x)\ > 2, or 

(I4b) To has no periodic critical points, except possibly that or 1 may 
be fixed points. 

Because To may be bi- valued at points in Co, a critical point c^o is considered 
periodic if there exists n > such that c^o € Tq^c^o). Condition (14) is 
necessary in order to ensure that the perturbed systems defined below satisfy 
uniform Lasota-Yorke estimates. Since we cannot exclude the possibility of 
the forward orbit of a critical point containing other critical points, these 
uniform estimates do not follow directly from the original paper [LY73] , but 
rather from later works, see §4.2. 

For what follows, we consider C 2 -small perturbations T e : / Q of To for 
e > 0. This means that a critical set for T e may be chosen as C t = {0 = 
co,e < ci,t < • • • < Cd : e = 1}, where for each i, e i — > Cj )£ is a C 2 function for 
e > 0. Furthermore, there exists 8 > such that for all sufficiently small e, 
there exists a C 2 extension T j£ : [c^o — 6, Cj+i 5 o + 5] — > IR of T e |[ c . ejCl+1 j, and 
T j(E — ► T^o in the C 2 topology. We also assume: 

1 A piecewise expanding map T : I O with C = {0 = Co < Ci < • • • < cj, = 1} as a critical 
set is weakly covering if there is some N such that for every i, u£L T* ([d, Ci+i]) = I 
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(PI) Unique ACIM. 

For e > 0, T t has only one ACIM fi e , with density <p e := d(i e /dx. 



(P2) Boundary condition. 

The boundary point does not move, and no holes are created near the 
boundary; precisely, 

(P2a) If b ^ Co, then necessarily To (6) = b. We assume further that for 
all e > 0, T € (b) = b. 

(P2b) If b G Co, we assume that To (&_) < b < T)(6+), and also that 
b £ C e for all e. 

If the boundary point does move under the perturbation, condition (P2) 
often can be satisfied by performing a small change of coordinates; see §2.4. 



2.2 Main results 

The central question of this study is, for small e, how can we asymptotically 
approximate /i e by a convex combination of /U/ and [i r l To that end, let 
H he := Ii n T~ l {I r ) and F r ,e := T D T" 1 ^). We refer to these sets as 
holes. Once a T e -orbit enters a hole, it leaves one of the invariant sets for To 
and continues in the other. As e — > 0, the holes converge (in the Hausdorff 
metric) to the infinitesimal holes from which they arise. 

Condition (PI) ensures that for e > 0, at least one of the holes has 
positive Lebesgue measure. In view of (13), without loss of generality, we 
suppose that m(Hi j6 ) > and define 

l.h.r. = lim 

if the limit exists, {l.h.r. stands for limiting hole ratio.) 

Theorem 1 (Approximation of the invariant density). Consider the 
family of perturbations T t of To under the assumptions stated in §2.1. Sup- 
pose that l.h.r. above exists. Then as e — > 0, 

L 1 OL 

4> e — > a<f>i + (1 — a)4> r , where = l.h.r.. 

1 — a 

We allow for l.h.r. = +oo, in which case a = 1. Several straightforward 
generalizations of the above result are discussed in §2.4. 

Remark 1. The limit l.h.r. will always exist as long as the perturbations 
open up holes Hi t whose size is truly first order in e: For simplicity, suppose 
that there are only two infinitesimal holes, hi £ Ii and h r G I r . Then we 
can always write e = (h* — a*e + o(e), h* + b*e + o(e)) for * G {/, r}, and 
if ai + bi > 0, then 

_ 4> r (h r )(a r + b r ) 
(f>i(hi)(ai + bi) ' 
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For example, this will be the case if T t = Tq + eg + o(e) for some smooth 
function g with g(hi) > 0. 

Remark 2. An alternative definition of l.h.r. is as a limit of a ratio of escape 
rates: For * G {/, r}, we can consider a dynamical system with a hole, where 
orbits stop upon entering the hole, by using the unperturbed map To|/ t with 
H* )fi as the hole. See [DY06] for an exposition of such systems. Let R* te be 
the exponential escape rate of Lebesgue measure and suppose that there is 
only one infinitesimal hole in each initially invariant interval. Then as e — > 0, 
/i*(iJ*, e )/i?*, e -» 1. [BY08, KL] 

Next, let C € be the Perron- Frobenius operator associated with T e acting 
on the Banach space BV = {/ : I — ► C : var(/) < oo} 2 with the variation 
norm, and let cr(C e ) denote the spectrum of C e . It follows from e.g. [Kel89, 
Thm. 8.3(b)] that Cq has one as an isolated eigenvalue of multiplicity two. 
Furthermore in [KL99] the authors show that for fixed small 5 > and for 
every e > small enough, cr(C e ) n Bg{l) consists of exactly two eigenvalues, 
1 and p e < 1, each of multiplicity 1. As e — > 0, p e — > 1 and the total spectral 
projection of C € associated with cr(C e ) n B$(l) converges (at a given rate 
in an appropriate norm) to the total spectral projection of Co associated 
with <t(£o) H B$(l). Note that in §4.2 we will show that the assumptions of 
[KL99] are satisfied in our context. 

Theorem 2 (Characterization of the eigenspace corresponding to 
the lesser eigenvalue). For each e > small enough, there is a unique real- 
valued function ip e S BV satisfying C e tp e = p e ip e , I^cIl 1 = 1; an d fj ipedx > 
0. As e -► 0, 

A — > -4>i - 

Remark 3. Suppose pi and p r are both mixing for To. Given a typical initial 
density / G 5V (i.e. one with nonzero coefficient of ip € when expressed as 
a linear combination of eigenvectors), as £™f — ► e , the deviation — 
e becomes approximately proportional to ip e for n large. In this case, 
Theorem 2 implies roughly that for e small, C™f becomes close to a linear 
combination of 4>i and <p r more quickly than it comes close to the specific 
linear combination afy + (1 — a)(j) r . 

2.3 Examples 

The three piecewise linear maps shown in Figure 3 satisfy assumptions (II)- 
(14) from §2.1. In all three cases, normalized Lebesgue measure restricted to 
the left or right intervals is the unique ACIM of the corresponding restricted 
system. 

2 Technically, two functions of bounded variation are considered equivalent if they differ 
on at most a countable set. Here and elsewhere, we generally ignore such distinctions. 
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hi b h r hi b h,. h, 

Figure 3: Piecewise linear maps giving rise to metastable systems. 



Adding a small C 2 perturbation g : I x [0, 6q) — > I such that g{-, 0) = 
and for e ^ 0, g(b, e) = 0, g(hi, e) > and g(h r , e) < gives a one-parameter 
family of perturbations T e := Tq + g(-,e) satisfying assumptions (PI) and 
(P2)- 

If lim e ^o Leb(H^ ,,! ) = l-h- r -i by Theorem 1, the invariant densities cj) e 
associated to T t satisfy 

L 1 Ct 

4> £ — > aLeb|/ ; + (1 — a)Leb|/ r , where = l.h.r.. 

The possibility l.h.r. = oo is allowed, and in this case, 

4> € — > Leb|/ r 

Other initial maps To for which Theorems 1 and 2 are applicable are 
shown in Figure 4. 





Figure 4: Examples of initial maps To which give rise to metastable systems 
for which our results hold. 



2.4 Generalizations 

Our results extend, with essentially the same proofs, to yield the following 
straightforward generalizations. 
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Multiple invariant sets. 



We can also allow To to have m > 2 invariant sets, provided it has a unique 
ACIM on each of them. The invariant sets may be intervals or a union of 
intervals. See Figure 5. In this case, the unique invariant density 4> t of T e 
converges to a convex combination of the initial ergodic invariant densities 
as e — > 0. The coefficients may be determined from m — 1 linear equa- 
tions involving limits of the quotients of measures of appropriate holes. If 
m > 2, the analogue of Theorem 2 says only that the eigenfunctions for 
T e whose eigenvalues approach 1 limit on the space of eigenfunctions for To 
with integral 0. 
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Figure 5: Initial maps To which give rise to metastable systems for which 
our results can be generalized. 



Boundary condition. 

The restriction that the boundary point does not move when To is perturbed 
is inessential; when it is relaxed, it simply means that the metastable states 
for T e are slight perturbations of the initial invariant sets. In this case, a 
smooth change of coordinates restores the hypothesis (P2). For example, 
when b ^ Co assumption (P2a) is actually superfluous, although the defini- 
tions in the statement of Theorem 1 must be modified slightly. As remarked 
earlier, necessarily To (6) = b. Furthermore, the graph of To intersects the 
diagonal transversely at this point. Thus for all small e > 0, there is a 
unique point b e near b satisfying T e {b t ) = b e . Then the quasi-invariant sets 
for T t are Jj j£ := [0, b e ] and := [6 e , 1], and the corresponding holes are 
defined by Hi )€ := Ii >e n T~ l (I r ^) and H r>e := I r , e n T" 1 ^). Aside from 
these minor modifications, the statements and proofs of our main results 
remain the same. 

When 6 G Co, (P2b) can be relaxed by no longer requiring that b £ C e 
for all e. In this case, when e > 0, C t contains a point b e that converges 
to b as e — ► 0, and the quasi-invariant sets and holes must be redefined 
as above. However, it is still essential to assume that no holes are created 
near the boundary, which we enforce with the assumption that Tq(6_) < 
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b < To(6 + ). For example, if T e (x) = [(3x mod 1/2) + 3e] • l x< \/2 + [(— 3x 
mod l/2) + l/2 — e] ■l x >i/2i then all of our assumptions aside from (P2) hold, 

with 6 = 1/2,//*= Lebl/,,, and l.h.r. = 1/3. However, as e — > 0, $> e — ► 
The difficulty is that orbits ejected from I r by T t immediately return to I r . 



Multiple limiting densities. 

When the limit l.h.r. in §2.2 does not exist, we let 

, , v ■ f PriHr e) — fJ, r (H r e ) 

l.h.r. = limmt — -, l.h.r. = limsup — r-. 

Since the function ^ r ^ r ' £ ? is continuous in e > 0, our arguments show that 
the set of limit points for <j> e as e — ► is precisely 

d 



a4>i + (1 — a)(j) r : G [ l.h.r. , l.h.r.} > . 

1 — a 



3 Proofs of the main theorems 

In this section, we briefly state the main properties of the invariant densities 
that will be needed before presenting the proofs of Theorems 1 and 2. For 
notational convenience, we will assume that there are only two infinitesimal 
holes, hi E // and h r £ I r ; the proof without this restriction is essentially 
unchanged. 



3.1 Properties of the invariant densities 

Here we record some of the relevant characteristics of the density functions 
(j) e ,(pi,(j) r . First, if / G BV, we can - and will - choose a representative of 
/ with only regular discontinuities, i.e. for each x, f(x) = (lim^^- f(y) + 
lim y ^ x + f(y))/2. Then, following [Bal07], we can uniquely decompose / = 
f reg + f sal into the sum of a regular and a singular (or saltus) part. Here 
jreg j g continuous with var(/ refil ) < var(/), and f sal is the sum of at most 
countably many step functions. We write f sal = J2ues s uH u , where S is the 
discontinuity set of /, s u is the jump of / at u, and H u (x) = — 1 if x < u, 
— \ if x = it and if x > u. This representation imposes the boundary 
condition f sal (l) = 0. Furthermore, vai(f sal ) = X^es \ Su \ — var (/)- 

Proposition 1 (Key facts about the invariant densities). There exists 
eo > such that: 

(i) Uniform bound on the variations of the invariant densities. 



sup var( 

0<e<e 



< +0O. 



Also, var(</>z), var(</v) < +oo. 
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(ii) Uniform bound on the Lipschitz constant of the regular parts. 

For < e < eo, each of the <p r e e9 is Lipschitz continuous with constant 



Also, <p r ^ 9 and <p r r eg are Lipschitz. 

(hi) Approximate continuity near the infinitesimal holes. 

For * G {l,r}, for each rj > 0, there exists 5 > such that for all 
< e < e , 

v&r[f lt -8,h t +8]{ < l )S e al ) '■= the variation of (p s e al over [/i* — 5, h* + 5] < r\. 
Also, <p* is continuous at /i*. 
The proof of Proposition 1 is technical, and so we defer it until §4.2. 

3.2 Proofs 

We recall that for any C X ,C 2 > 0, {/ G BV : \f\ Ll < d,vax(f) < C 2 } is 
pre-compact in L 1 . This fact will be used repeatedly in what follows. 

Proof of Theorem 1 

Using (i) of Proposition 1, we are able to choose a sequence of values e' 
converging to such that (p e > converges in L 1 to some function, which we 
denote by (pQ. Using the fact that (p e is a fixed point of the Perron- Frobenius 
operator C e associated to T e (see §4.1 for the definition), one can verify 
that 0o is an invariant density for To, and so there exists a such that (po = 
a(pi + (1 — a)(p r . We will verify that necessarily a/(l — a) = l.h.r.. From this 
it follows that there is exactly one limit point of <p e as e — ► 0, and Theorem 1 
follows. 

Now for e' > 0, fi e i(Hi e i) = ij, e i(H re /), because (p e / = d/j, e </dx is an 
invariant density for T e i. We will show that as e' — * 0, 



from which the equation a/(l — a) = l.h.r. and hence Theorem 1 follows 
immediately 

We prove only Equation (1), since the proof of Equation (2) is analogous. 
Write 



Lip((p T e eg ), and 



sup Lip((p r e e9 ) < +oo. 



0<e<e 



Ve'{Hl,e>) 
V-e'{Hr,t') 



OlfJ-l(Hi^) + o(l) • l2i(H lje <), 

(1 - a)^ r (H r ^) + o(l) • pL r (H r ^) 



(1) 

(2) 
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But as e' — > 0, H[ e i — > hi in the Hausdorff metric, and then /j,i(Hi te i) /Leb(Hi >e i) 
<f>i(hi) > 0, because cfri is continuous at hi. Thus our proof is completed by 
the following: 

Lemma 1. As e' — > 0, 

sup \(f> e i(x)-at<f>i(x)\—>0. 

Although this uniform convergence might at first seem surprising, Propo- 
sition 1 (ii) and (hi) essentially say that near hi, {(j> e >} behaves like a family 
of equicontinuous functions. 

Proof. We proceed by contradiction. Suppose that there exists C > and a 
subsequence e" — > of the e' values such as that for each e" , there is a point 
x e " £ Hi^n with \(f> e ti{x e ii) — ai4>i(x e ")\ > C. Necessarily, x e » — ► /i; as e" — ► 0. 

We restrict all functions of interest to the left subinterval Ii. Set j e n := 
4> r J, g — a(fi e9 and w e » := 4> s t % 1 — acf)f al , so that e » —afy = j € "+oj e 'i. Using (ii) 
of Proposition 1, let L be such that for all sufficiently small e" , Lip(7 £ ») < L. 
Next, we use (iii) with rj = C/5 and make sure to choose the corresponding 
5 < C/(5L) small enough so that \SLi^ hl _^ hl+ ^{a(f)f al ) < C/5 as well. Thus 
v&Ylhi-S^i+sfae") < 2C/5. Then if x € [/t; — <5, /i; + 5] , and e" is sufficiently 
small, x e » S [hi — 5,hi + 5] and 

|7 e »(x) +U e »(x)\ 

> |7 e "(x e ») + ^e" i x e") I ~~ l7e"( x ) +W e »(x) — 7 e ''(x e «) — a; e //(x e //)| 

> C - [L • 2<5 + 2C/5] > C/5. 

But this contradicts that 7 e « + uv = 4> e /> — a<fii — ► 0. □ 
Proof of Theorem 2 

First, we observe that the results of [KL99] guarantee that for small e > 0, 
p e < 1 is a simple eigenvalue of multiplicity 1 . Hence there are exactly two 
real-valued eigenfunctions, ±ip t , satisfying C e tp e = p e ip f and (V^Il 1 = 1- But 
for such functions, J ip e dx = J C e i() e dx = p e J ip e dx, so f ijj e dx = 0. We 
have the following uniform bound on their variations, whose proof we defer 
until §4.2. 

Lemma 2 (Uniform bound on the variations of the ip e ). There exists 
t\ > such that supo <e < ei v&i(ip t ) < +oo. 

Let be any limit point in L 1 of ip e as e — ► 0. Then, since p e — ► 1, 
it follows that ipo is invariant under Co, and is thus a linear combination 
of <pi and 4> r . Since l^ol^i = 1 and fiftodx = 0, necessarily ipo = ±|<^ =F 
\<t> r - Hence we can uniquely specify ip e by the condition fj ip e dx > 0, and 
Theorem 2 follows. 



13 



4 Proofs of the properties of the densities 



In order to prepare for the proofs of Proposition 1 and Lemma 2, it will 
be convenient to first show how to derive such properties for an invariant 
density of a single, fixed piecewise expanding map. We do this in §4.1. 
Then, in §4.2, we prove Proposition 1 and Lemma 2 by showing how such 
estimates can be made uniformly for the family of maps T e , e > 0. 

Before beginning, we remark that if / G BV, then for each x, |/| < 
\f(x)\ +var(|/|) < \f{x)\ +var(/). Integrating, we find that (/(^ < \f\ L1 + 
var(/). We will use this fact repeatedly below. 

4.1 Properties of an invariant density for a single piecewise 
expanding map 

Let T : I Q be a piecewise C 2 uniformly expanding map, with C = {0 = 
co < c\ < ■ ■ • < Cd = 1} as a critical set. Let £ be the associated Perron- 
Frobenius operator, i.e., the transfer operator acting on densities. We begin 
by briefly reviewing a method for finding an invariant density of T. Such a 
method was introduced in [LY73]; see Chapter 3 in [BalOO] for a more modern 
exposition. Let Ay = inf^nc" l^'C^)! > 1 be the minimum expansion and 
D T = swp x <£ C \T"(x) \ / \T'(x)\ be the distortion of T. Then if / G BV, 
x $ TC, 

d 

Cf{x) = Y,f&{xWi{x)\^J i (x), (3) 

i=l 

where J { = T| [c ._ 1)C .]([ci_i, Cj]) and & = (Tl^^^) -1 : Ji -> One 
can show that there exists constants (3 G (0, 1) and Cly such that for each 
n > 1 and / G BV, the following Lasota-Yorke inequality holds: 

var(£ n /) < C LY /3 n var(/) + C LY \f\ L1 . (4) 

In fact, P can be chosen as any number greater than A^ 1 , although we will 

not use this fact. Set F n = ± Y.l=l Then F n ^ where 4> G BV 
is the density of an ACIM for T. Using Helly's Theorem, one has that 
var(<£) < C L y- 

We wish to characterize the properties of the regular and singular terms 
in the decomposition <fi = <p re9 + cf) sal . First, let us define a hierarchy on the 
set of points in the postcritical orbits S = Uk>iT k C by := m£{k > 1 : 
u G T k C}. The following characterization is motivated by the discussion of 
the invariant densities for unimodal expanding maps found in [Bal07] and 
[BS08]. In particular, in [BS08, §3.3] a norm is introduced on the sequence of 
jumps of 4> along the postcritical orbit with weights that grow exponentially 
in 

Lemma 3. Given the hypotheses above, 
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(a) <ff e9 is Lipschitz continuous. Furthermore, there exists a constant 
C dis = C Ais (\ T ,D T ) such that Lip(^) < C dis (l+C LY ). C dis (^T,D T ) 
can be defined so that it depends continuously on At > 1, Dt > 0. 

(b) The discontinuity set of ^ is a subset of S = U k>1 T k C. If we write 
sal = J2ueS s uH u , then for each m > 0, E{ uei S:#(«)>m} W < 

\ — m s~t 
A T OlY- 

Proof. We begin by noting from Equation (3) that F n is smooth except 
possibly at points in U^ _1 T fc C. Write F° al = Y^ u eS s u, n H u . Then we can 
show that \s u n \ decays uniformly exponentially fast in i.e. 

Sublemma. For each m, n > 0, E{u&S: #(«)>m} K,n| < A^ m C L Y- 

Proof o/ i/ie Sublemma. If m > n, ^#(n)>m l s «,",| = 0, and Xl#(«)>o I s ",™! = 
var(F,f ') < var(F n ) < C LY . Since F n = ^£F n -l + if #(«) > 1 we see 
from Equation (3) with / = F n _i that \s %n \ < ^A^ 1 1]{ 1 , eiS:T , ;=u } K,n-l|- 
Thus if < m < n, 

^ l S «:"l — X/ « ^T" 1 1 1 

#(u)>m #(u)>m {ue5:Tt)=u} 

< A T > Su,n— 1 S " " " 



< 



#(u)>m-l 

n — m 



^t" 1 ^ ] | s u,n— m| — Aj. m CLY- 



n 

#(«)>0 

In the inequalities above, we use the fact that if #(u) > 1, then T~ l (u) does 
not contain any critical points. □ 

Using a diagonalization argument, we may find a subsequence rij such 
that for each u, s u ^ nj converges as rij — > oo to some number, which we write 

as s u . In particular, for each m, Xy#(u)>m 

| < A^ m C LY , and F° al 
psai ^ w j iere we define F sal = ^2 u& s s u H u . Furthermore, a standard distor- 
tion estimate (see for example the proof of Proposition 3.3 in [Bal07]) shows 
that there exists a constant Cdi s such that for each n > 0, Lip((£ n l) re9 ) < 
Cdis |£ n l|oo — Cdis(l + var(£ n l)). Here, C^m depends only on the minimum 
expansion and on the distortion of T. In particular, sup n>1 
Cdis(l + Cly)- By the Arzela-Ascoli Theorem, we may find a continuous 
function F reg such that some subsequence of {F^ 9 } converges in L°° to 

preg _ 

By the uniqueness of the decomposition (f> = (ff e9 + cf> sal , we conclude 
that (f) re 3 = F re 9 and (j) sal = F sal . Lemma 3 follows. 

□ 
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4.2 Proofs of Proposition 1 and Lemma 2 

We prove only the claims about (J) e for e > 0, and leave the claims about 
(pi, </y to the reader. 

Let C e be the Perron- Frobenius operator (3) associated to T e . The first 
key step is to prove that the C e with e sufficiently small satisfy Lasota- 
Yorke inequalities with uniform constants. Let A e and D e be the minimum 
expansion and distortion of T e , respectively. Then as e — > 0, A e —* Ao and 
D e —* Z?o- Furthermore, T e is a piecewise C 2 uniformly expanding map 
that is a small C 2 perturbation of To, and the two critical sets C e , Co are 
e— close together. This is not sufficient to guarantee uniform Lasota-Yorke 
inequalities, see for example [Kel82, §6] or [Bla92]. However, such uniform 
inequalities do follow with the additional assumption (14), which guarantees 
that either (a) we have Ao > 2 or (b) Tq has no periodic critical points, 
except possibly the points in dl as fixed points. We assume the former case 
in our presentation here, and comment on the latter case at the end of this 
section. 

Fix A £ (2, Ao). The original proof from [LY73] shows that if / G BV is 
real- valued, 

var(£J)<(2A- 1 )var(/) + C e |/| L1 , 

where 

C e = D e /\ e + 2 max \q+i >e - Ci J" 1 . (5) 

i 

(Compare also [Liv95][§2].) Iterating, we find that for sufficiently small e, 
for all such / and n > 1, 

var(£™/) < /Tvar(/) + C LY \f\ Ll , (6) 

with/? = 2A _1 and Cly = 2Cb/(l— 2A _1 ). Similar estimates can be made for 
complex-valued / by applying (6) to the real and imaginary parts separately. 
Since each T e has a unique ACIM, we know from our discussion in §4.1 that 

for sufficiently small e > 0, i ^Cfc=o ^ ^e asn ~ > 00 • It follows from 
Lemma 3 that var(</> e ) and Lip((/£ e9 ) are uniformly bounded. 

Next, we prove (iii). Given ij > 0, choose n large enough that A~ n Cx,Y < 
rj. Using (12), we can choose 5 > so small that for < k < n, (TqCo) n 
[K - 25, K + 26} = 0. It follows that for e sufficiently small, (T e fc C £ ) n [h* - 
S, h* + S) =0 as well. Using part (b) of Lemma 3 with m = n, we then see 

that yaT[ ht -s,h t +8](€ al ) < V- 

Finally, to prove Lemma 2, we use Equation (6) with f = tfj e , n chosen 
so large that (3 n < 1/2, and e chosen so small that p" > 3/4. It follows that 
var(VO < C LY M ~ P") < ^ly- 
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Modifications when the minimum expansion is not bigger than 
two 

If, in assumption (14), the minimum expansion of To is Ao < 2, one derives 
Lasota-Yorke estimates for Cq by first fixing TV large enough so that > 
2. Then the arguments from [LY73] used above will yield a Lasota-Yorke 
estimate for Cq , and this can be interpolated to give similar estimates for 
Cq. One can try to obtain uniform estimates for C e , but the arguments 
used above will only work if the critical points for are in a one-to-one 
correspondence with and very close to those of ', compare Equation (5), 
as would be the case if C n (U^T k C ) = 0. 

7XJMZ 

To If T e if 

Figure 6: Creation of small intervals of differentiability. 




Unfortunately, this will never be the case in our setting, at least when 
b E Co- This is because the infinitesimal holes in Hq are necessarily critical 
points, and they are mapped to b by To. Because at least some of the 
infinitesimal holes must be mapped across the boundary point when e > 0, 
this means that necessarily T e 2 will have more critical points than Tq, and 
these additional critical points will create very short intervals on which T e 2 
is smooth; see Figure 6. However, this problem can be dealt with using 
assumption (I4b). Specifically, in [BY93] it is shown that because of the 
restriction on the periodic critical points, the growth in the number of the 
very short intervals on which T e n is smooth as n increases can be controlled, 
and that uniform Lasota-Yorke estimates can still be made. Precisely, there 
exists eo > and constants (3 E (0,1), Cly such that for each e E [0, eo], 
n > 1 and / E BV, 

var(£-/) < C LY /3"var(/) + C LY \f\ Ll . 

The proof of this is essentially identical to the proof of Lemma 3.2 in [BalOO] 
(see also her Remark 3.4, and compare the proof of Lemma 8 in [BY93]), 
and so we omit it. The rest of the proofs of Proposition 1 and Lemma 2 
proceed as above. 
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